Spectral shape analysis of ultraviolet luminescence in n-type ZnO:Ga 
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Thin films of laser molecular-beam epitaxy grown n-type Ga-doped ZnO were investigated with respect to 
their optical properties. Intense room-temperature photoluminescence (PL) in the near-band edge (NBE) region 
was observed. Moreover, its broadening of PL band was significantly larger than predicted by theoretical results 
modeled in terms of potential fluctuations caused by the random distribution of donor impurities. In addition, 
the lineshape was rather asymmetrical. To explain these features of the NBE bands, a vibronic model was 
developed accounting for contributions from a series of phonon replicas. 

PACS numbers: 78.55.Et, 81.15.Fg, 71.35.Cc, 72.15.-v 



Recently, optical properties of near-band-edge (NBE) 
recombination in ZnO epitaxial layers have been extensively 
investigated. Compared to undoped ZnO, the properties of 
heavily-donor doped ZnO, such as phonon interactions, are 
not well understood. In our previous work, we reported 
the observation of intense NBE photoluminscence (PL) for 
ZnO:Ga at room temperature (RT). Figure 1(a) shows their 
PL spectra at four different doping levels (dashed lines). 
The experimental procedures yj] were identical to those 
adopted in our previous study and their details were reported 
elsewhere |'2]. The shape of the PL band is asymmetrical. 
The linewidths of the PL increase from 154 to 195 meV 
with an increase in Ga concentration {nca)- Previously, we 
estimated the "theoretical" broadening for the PL in terms 
of potential fluctuations caused by the random distribution 
of donor impurities (cf. Fig. 3 of Ref. 1). These values are 
also listed in Table I. However, the calculated widths (cti) 
were smaller than the experimentally observed broadening. 
Table I: Sample specification and main characteristics for the 
four n-ZnO samples. Ga concentration (nca), theoretical 
broadening (Ti, broadening evaluated from the Monte Carlo 
simulation a-2, zero-phonon peak energy E, and the deduced 
Huang-Rhys factor S. 
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I'hese theoretical widths might be exceedingly underes- 
timated because the model might be too much simplified. 
We here adopt another approach properly taking the micro- 
scopic fluctuation of donor concentration into account; the 
thermalization redistribution model developped by Zimmer- 
mann etal\3]. According to this model, there is a relationship 
between the Stokes shift and the full-width at half maximum 
(FWHM) of PL. One can estimate the energy scale of the band 
potential profile fluctuation from the Stokes shift and predict 
consistently the FWHM. Figure 2 shows a Stokes shift of the 
luminescence plotted against the gallium concentration. Ob- 
viously, the shift energy increases with incorporation of the 
donor dopants. 

At sufficiently high temperatures, the above-mentioned re- 
lationship is indicated as Estokes = —cr^/kT. For exam- 
ple, the Stokes shift for the sample with the nca of 8.0 
x 10^^ cm^'^ is 22 meV, leading to a of 24 meV. It also yields 
the temperature independent value of 2cry^ln4) ~ 46 meV 
in the FWHM. This can be also supported by the Monte Carlo 
simulation. We now simulate a PL spectrum. 

Miller- Abraham's rate for phonon-assisted exciton tunnel- 
ing between the initial and final states i and j with the energies 
of Ei and Ej was adopted to simulate the hopping events of 



excitons or of photocreated carriers: 



vq exp 
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Here rij is the distance between the localized states, a is the 
decay length of the wave function, and vq is the attempt- 
escape frequency. Hopping was simulated over a randomly 
generated set of localized states with the sheet density of N. 
Dispersion of the localization energies was assumed to be in 
accordance with a Gaussian distribution. 
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FIG. 1: (Color online) (a) Room-temperature photoluminescence 
spectra (dashed curves) of n-type ZnO doped with different Ga con- 
centrations. Also shown are the results of the fit to the data using a 
'vibronic' model (solid lines), (b) The lowermost PL curve of frame 
(a) with individual contribution of the emission lines (dash-dotted 
lines). 
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FIG. 2: The PL Stokes shifts (closed squares) plotted against the Ga 
concentration. 
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FIG. 3: Simulated PL spectrum under stationary excitation. The Ga 
concentration is shown in the figure. 
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With the peak positioned at the mean excitonic energy Eq and 
the dispersion parameter (the energy scale of the band poten- 
tial profile fluctuation) a. For each generated exciton, the hop- 
ping process terminates by recombination with the probability 
Tq" ^ and the energy of the localized state, where the recombi- 
nation has taken place is scored to the emission spectrum. 

In addition, the thermal broadening effect was implemented 
by convoluting the above spectrum with a Gaussian curve, 
which describes the thermal part with a width of LSfc^T \^. 
The nomenclatures (ks and T) take their conventional mean- 
ings. Because the above equation contains the term of ther- 
mal activation, at a glance, it seems that it is unnecessary to 
perform the convolution. However, because Kazlauskas et al. 
adopted this treatment, we follow it in this work J^. Mean- 
while, the overall FWHM of the PL band is given: 

FWHM (aa) ~ [i2ay/{lnA))^ + {1.8kBTf]^^\ (3) 

The solid line in Fig. 3 represents the simulated result ob- 
tained for the following values of parameters: Na^ — 0.1, 
tqI'o = 10^, and cr ~ 24 meV. The line-shape of resulting 
spectrum was independent of the choices of the Na^ and tqI'o 
parameters. The simulated curve could be well approximated 
by the Gaussian function, which is not surprising if consider- 
ing the measurement temperature (300 K). We summarize the 
FWHMs ((T2) of the simulated PL spectra for the samples at 
four different doping levels in Table I. 

The simulated broadening ((T2) is still smaller than that of 
experiments. This discrepancy is thought to be due to the con- 
tribution of the phonon replicas. In addition, the simulated 
line-shape is rather symmetrical, which is in great contrast 
to the asymmetrical shape observed in the experimental spec- 
tra. In compensated semiconductors, even a low-temperature 
emission spectrum shows very broad feature. Its line-shape is 
sometimes asymmetrical. Kuskovsky and his co-workers | ^ 
have submitted a theory reproducing such an asymmetrical 
line-shape. On the other hand, our experiments were per- 
formed at room temperature. Therefore, it is natural to think 
that the shape of a unresolved zero-phonon band is rather sym- 
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metrical because of the thermalization effects. We did not 
estimate the width and the asymmetry of our PL using their 
theory because, unfortunately, it is applicable only to the case 
of extremely low temperatures and because our samples are 
not very severely compensated We think that the zero- 
phonon peak and its phonon replicas are not spectrally re- 
solved well in the spectrum because intensity of a one-phonon 
replica is comparable to that of a zero-phonon band, leading 
to the larger experimental broadening. 

The typical consequence of the interaction with phonons is 
the appearance of phonon-assisted emissions in a PL spec- 
trum. We infer that the coupling with longitudinal-optical 
(LO) phonon is very efficient in this highly polar material. 
The reported coupling constant of ZnO is nearly four times 
stronger than that of GaN and about eight times stronger than 
thatofZnSe |7|. 

For simplicity, we took only the terms of longitudinal opti- 
cal (LO) phonons into account 1'8','9','To], whose interaction is 
the strongest. In this case, the transition energy of the repli- 
cas will be hu = EzpL — Vt^oE^O, where i?LO=72 meV in 
ZnO, Ezph is the energy of the zero-phonon peak and ?7lo is 
an integer. 

Each individual emission line in the spectrum was modeled 
by a Gaussian function with linewidth ai (cf. Table I), which 
can be justified by the Gaussian-like lineshape observed in our 
simulated spectrum. The values of <ti were determined from 
the above-mentioned theory because these values are larger 
than those of the Monte Carlo simulation, 0-2. The probability 
of a given phonon emission is proportional to (S'lo'^/^/lo!)' 
where 5lo is the Huang-Rhys factor for the LO phonons and 
r^LO the number of phonons emitted in a transition. The solid 
lines in Figs, ^a) and (b) were obtained by summing all of 
the emission lines using an appropriate set of parameters. The 
donor-impurity concentrations (nca), zero-phonon peak en- 
ergies (E), and S'-factors are compiled in Table I. Five dash- 
dotted curves in Fig.^^b) denote the respective contributions 
whose peak positions are equidistant by the E^lo of 72 meV. 
The important conclusion from this fit is that the asymmetri- 
cal and broad PL band can be explained by the contribution of 



LO phonon replicas. It should be noted that the Huang-Rhys 
factor (S'lo) depends on the ionic nature of semiconductors 
and on the separation of Fourier transformed charge distribu- 
tions. This sometimes deviates from the bulk value [8, 9,131. 
Very recently, Shan et al. 1 1 1] reported that the RT PL spec- 
trum in undoped ZnO is dominated by the phonon replicas of 
the free exciton transition with the maximum at the first LO 
phonon replica. This corresponds to the S'-factor exceeding 
unity. In addition, in the undoped ZnO, the zero-phonon emis- 
sion becomes weaker than ILO-phonon-assisted annihilation 
processes at elevated sample temperatures. 

Even if the values of ai are comparable with the optical 
phonon energy E^o^ due to the high measurement temper- 
ature, we believe it is probable that the radiative transition 
would go via phonon-assisted path while the impurity levels 
are available at the same energies. 

We see the shoulders at the ILO phonon replica position 
in both the fitting and experimental results, particularly in the 
first two cases. The "experimental" shoulder is less prominent 
than that of fit. The reason of this is not clear at this moment 
but we speculate that this is because of our assumption that 
the width of the phonon replicas is same with that of the zero- 
phonon band. We think that the broader replicas are more 
realistic. 

In summary, gallium was successfully applied as a source 
of donor impurities in the laser molecular-beam epitaxy 
growth of n-type ZnO. We presented spectroscopic measure- 
ments to study the evolution of spectra in «-type ZnO along 
with Ga concentrations. Asymmetric and broad photolumi- 
nescence spectra were fitted to a model that takes into account 
phonon replicas, in which the linewidths of each individual 
emissions have been determined mainly from the concentra- 
tion fluctuation of dopants. 
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